Heisenberg Spins on a Cylinder Section * 
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Classical Heisenberg spins in the continuum limit (i.e. the nonlinear cr-model) are studied on an 
elastic cylinder section with homogeneous boundary conditions. The latter may serve as a physical 
realization of magnetically coated microtubules and cylindrical membranes. The corresponding rigid 
cylinder model exhibits topological soliton configurations with geometrical frustration due to the 
finite length of the cylinder section. Assuming small and smooth deformations allows to find shapes 
of the elastic support by relaxing the rigidity constraint: an inhomogeneous Lame equation arises. 
Finally, this leads to a novel geometric effect: a global shrinking of the cylinder section with swellings. 
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I. INTRODUCTION 

Microtubules^] and cylindrical membranes abound in 
nature. They are important both in biological context 
and in industrial applications. The coupling of mag- 
netism with the deformability of these soft materials 
opens up new, exciting avenues of investigation. We can 
envision these microtubules to be synthesized from cither 
magnetic organic materials or poiymersa or to enclose 
various magnetorheological nuidscl in such a way that 
the magnetic and elastic properties are coupled. Specifi- 
cally, the shape deformation can be interpreted in terms 
of bending elasticity of the membrane. To explore mag- 
netoelastic properties we will approximate the surfaces of 
microtubules and cylindrical membranes as a continuum 
of classical Heisenberg-coupled spins. 

Our motivation is to amplify on the concept of ge- 
ometrical frustration, i.e. mismatch of different length 
scales in a system leading to deformation of the geome- 
try. The length scales of interest are the characteristic 
width of the spin (magnetic) soliton and the underlying 
geometric length of the cylinder (i.e., radius). The main 
idea is that the mismatch of length scales is reflected in 
the violation of the so called self-duality equations and 
the increase in magnetic energy of the system from that 
of the minimum energy in the corresponding homotopy 
class. Restoring self-duality or the minimum energy is 
achieved by deforming the manifold in such a way that 
geometrical frustration is relieved. 

The continuum limit of the Heisenberg Hamiltonian for 
classical ferromagnets or antiferromagnets, for isotropic 
spin-spin coupling is the nonlinear cr-modelETEl. The total 
Hamiltonian for a deformablc, magnetoelastically cou- 



where H m agn, 7~L e i and Hm-ei represent the magnetic, 
elastic and magnetoelastic energy, respectively. In the 
present paper we will focus on the magnetic part and 
the elastic part only since, for quasi-one-dimensional spin 
solutions on the cylinder, H. m - e i merely renormalizes 
7~L m agn- For the nonlinear cr-model, the magnetic energy 
onpa-curved surface S, in curvilinear coordinates, is given 
byiS 
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where J denotes the coupling energy between neighbour- 
ing spins. The order parameter n is the local magneti- 
zation unit vector specified by a point on the sphere S 2 . 
The metric tensors (gij) and (h a p) describe respectively 
the support surface S and the order parameter manifold: 
as customary, dO represents the surface area and g the 
determinant det(gij). 



II. RIGID CYLINDER SECTION: 
NON-SELF DUAL SOLITONS 

First, let us consider the nonlinear cr-model on a rigid 
cylinder section. 



A. Formulation 

For our purposes a suitable representation, in cylindri- 
cal coordinates (p, p, z), is given by 



pled manifold is given by Ti = Ti. r: 



Hel + TL 



m — eh 



(2) 



*Internat. J. Modern Phys. B 14, 2093 (2000) [cond-mat/9909095] 



1 



where the constant parameter r is real and positive, while 
the angle tp varies from — n to n. One can easily check 
that the metric is given by 

g = r 2 dp®dip + dz<g)dz, (3) 

therefore g tpz = g z,p — and we have 

9 VV V9=K 9 zz V9 = r. (4) 

From now on, we restrict ourselves to a section of the 
cylinder: the z variable will vary from — Az to Az where 
< Az < oo. 

As usual, the local magnetization n is described by 
its spherical coordinates (0,$), then the metric on the 
Heisenberg sphere is given by 

h = d9<g>d9 + sin 2 9 d$<g>d$. (5) 

Assuming homogeneous boundary conditions at both 
ends (9 = [ir] as z — > ±Az) allows to map each bound- 
ary of the section to a point: thus we compactify our 
cylinder section into the sphere S 2 . Consequently, the 
mapping of our support to the order parameter mani- 
fold is classified by the homotopy group II2 (S 2 ) which 
is isomorphic to Z: spin configurations may be classified 
according to their homotopy classD'Eil. 



where g = q v £. 

Equation @) may be integrated once to yield 

9 2 =sin 2 9 + m me[0,+oo). (10) 

Performing the change of variable sin.©| = dn(u | 1 + fh) 
where dn is a Jacobi elliptic functior£3, the differential 
equation becomes u 2 g = 1. Let us denote by a(- | fh) the 
increasing solution of ( pb| ) specified by the parameter fh 
and subject to the boundary condition a(0 | fh) = 
Readily, we get: 

sma(g \ fh) = dn(g \ 1 + fh) . (11) 

Thus, the general SG solution, in natural coordinates, is 

9(g) = ea{g \ fh) + c s = ±1. (12) 

Further the function a(- | fh) satisfies 

a(g + 2qK a \ fh) — a(g \ fh) + qir q G Z, (13) 



B. Derivation of non-self dual solitons 

Henceforth, without loss of generality, only cylindrical 
symmetric configurations (d v Q — d z Q = 0) will be con- 
sidered. Thus the magnetic Hamiltonian (|l|) becomes 



+ AZ +7T 

J J dzjdf 

-Az -ir 
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where a subscript stands for differentiation. Rescaling 
the z variable in equation (Jq) gives: 



Hmagn = J J d( J dip [6 2 
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where C = z/r and = Az/r. The Euler-Lagrange 
equations corresponding to (Q) are: 



$ w = 0, 

e^c = $ 2 sine cos e. 

From it follows that 

$ v = q v q v G Z. 



(8a) 
(8b) 

(9a) 



Substituting this into ( ph| ) and rescaling again the rotat- 
ing angle, we get the sine-Gordon (sg) equation 



<d ge = sin 9 cos 0, 



(9b) 



where the quasi quarter-period K a is relatecLto the com- 
plete elliptic integral of the first kind K bycj 



K a (m) = K(l+m). 



l =i4 1 



l+m 



(14a) 
(14b) 



Clearly, the parameter m tunes the quasi quarter-period 
K a . 

Using the solutions of equation (pb|), the g^7r-soliton 
configuration consistent with the boundary conditions 
can be obtained easily. Up to an irrelevant additive mul- 
tiple of 7T, we have 



9(C) = ea(q v ( \ m) - 5 e ven,q c f , 



where the parameter m is given by 



»> = K-A ^AC 



(15a) 



(15b) 



If we consider an infinite cylinder, the solution ([15]) rep- 
resents a soliton lattice with a period 2K a . When fh 
tends to + , this period tends to 00 and we recover the 
self-dual equations!^. 
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C. Geometric frustration 



only the elastic energy which is essentially stored in the 
bending of the deformable supportliTlia: 



The magnetic energy E maff „ of the above configuration 
( |l5| ) may be compared with the corresponding topolog- 
ical minimum energy E which does not depend on the 
geometry of the support-manifold. Performing the Bogo- 
mol'nyi's decomposition^ yields 



^gdQ(H- H ) 



(20) 



E = 8ttJ| 



(16) 



where the topological charge (i.e. the winding number) 
Q equals to qc,q v - Let £ m agn denote the ratio E mas „/E. 
A straightforward calculation shows that £ m agn depends 
on the parameter rh only; we have 



Here the constant k c denotes the bending rigidity of the 
(cylinder) material, H represents the mean curvatureE3, 
and the spontaneous mean curvature Hq tends to bias 
the mean curvature for recovering the spontaneous shape. 
Assuming small and smooth deformations and expanding 
to second order in A, A, and A,, lead to 
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(17a) 
(17b) 
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(21) 



which increases strictly from 1 to oo with respect to m; 
E is the complete ellipti c int egral of the second kindt3. 
Therefore, according to (15b) E ma gn is increasing func- 
tion of the radius r and decreasing function of Az. 



III. ELASTIC CYLINDER SECTION: 
FRUSTRATION RELEASE 

Next, let us relax the rigidity constraint and consider 
the nontrivial spin configuration on an elastic cylinder 
section. Accordingly, the soliton will try to minimize its 
magnetic energy E magn to th^Hainimum energy E by de- 
forming the elastic supportclil3't2l. 



A. Deformable metric 

Since the radius r appears as the relevant geometric 
parameter and z variable as the relevant curvilinear co- 
ordinate, we relax r with respect to z and write 



r(z)=r [1 + A(z)] 



(18) 



where rg represents the spontaneous radius and the func- 
tion A describes local deformations. The metric on the 
deformable manifold in this case remains orthogonal, and 
(El) reads 



g zz Vg= 



(19) 



B. Derivation of small and smooth deformations 

As the problem is quasi-one-dimensional, the magne- 
toelastic energy 7i m _ e ; merely renormalizes the-spin cou- 
pling energy J in the magnetic energy 7i mag „ll3. There- 
fore, we add to the nonlinear er-model Hamiltonian ([j]) 



Before deriving the Euler-Lagrange equation for the 
total Hamiltonian H = H m agn + 'H&i , we calculate the 
magnetic energy associated with the nontrivial spin con- 
figuration (|l5|). Expanding to second order in A and A z 
the relations (19) enable to rewrite (0) as follows 



n 



+ AC +7T 

magn= J J d( j dip 
-AC -w 



[e 2 + sin 2 e $ 2 ] 



(-A+1A 2 ) [9 2 -sin 2 e$ 2 ] 
+ A 2 sin 2 8 $ 2 



(22) 



On the other hand, according to ( p^| ) and (|IT|), the spin 
configuration (jlj) verifies 







c 



sin 2 6 $1 



(1 + m) q% [1 - sn 2 (^C I 1 + to)] 



\ — V- ' "V lip 

Inserting these relations into (|22| ) and simplifying, we ob- 
tain 



+ac 
2-Kjq% J dC 

-AC 



(2 + to) - to [-A 



(1 



[2 + A 2 ] sn 2 (^C | 1 + m) + A 2 



(23) 



The Euler-Lagrange equation for the problem takes the 
following form 

A eg + [(1 + m)A - mB sn\g | to)] A 

where we have set 

1 + to, 
2 L 1 c 2 " 



to 
.4 



B = 



mqfp 
2 

rag 2 ' 



1 + — 
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qlVT 



rn 



to j, (24a) 

(24b) 
(24c) 

(24d) 

(24e) 
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Here we have introduced the relative coupling energy 
k = J/k c and the variable g = q v Q (see (|9b|)). The 
l inea r inhomogeneous second-order differential equation 
(24a) is related to the well-known homogeneous (Jaco- 
bian) Lame's equation which occurs in several ph ysica l 
contextscjO. We have found no direct treatment of (24a) 
in the literature. However, an, approach based on the 
derivation of the Lame functionsEH allows to find a par- 
ticular solution denoted by L(- | m;A,B,j). Therefore, 
if small and smooth deformations are assumed, a suitable 
deformation function A is given by 



A(C) = L(<^ CI 1 + A, B,j) 



C. Frustration release 



(25) 



The geometric frustration becom es e vident in consid- 
ering the Euler-Lagrange equation (24a): A = (i.e. the 
rigid cylinder section) is not a solution. Now let us turn 
our attention to the deformation mechanism. As we have 
seen, the magnetic energy F, magn of the configuration (|l5| ) 
is an increasing function of the radius r: therefore for 
sufficiently small and smooth deformations, the soliton 
tries to decrease the radius r. In other words, the soliton 
tends to collapse the cylinder section. On the other hand, 
the elastic Hamiltonian ( p0| ) tries to maintain the spon- 
taneous shape of the deformable support. Consequently, 
the competition between the magnetic energy and the 
clastic energy induces a decrease of the radius r. Since 
according to @) the soliton energy is essentially localized 
in the spread zone, the deformation is less important in 
the spin flip region. 

To understand the geometric meaning of the frustra- 
tion release, let us define the relative dilatation A by 
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FIG. 1. The relative dilatation A as denned in (^6|) cor- 
responding to the deformation function (E5fi associated with 
cylinder sections in presence of a 27r-soliton versus C/^C for 
different values of AC and with the relative coupling energy 
k, fixed to 1/16. 
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(26) 



As shown in Fig. |l|, the decrease of the radius leads to 
a global shrinking whereas a local swelling arises where 
the spins twist. 



IV. CONCLUSION 

In conclusion, we showed that cylindrical magnetic 
surfaces, if approximated by classical Heisenberg spins, 
exhibit novel magnetoelastic deformation effects in the 
presence of spin solitons due to a mismatch of length 
scales. In particular, we showed that the Euler-Lagrange 
equation for the nonlinear cr-model on a cylinder sec- 
tion with homogeneous boundary conditions is the sine- 
Gordon equation. As announced, the frustration release 
leads to a novel geometric effect: the cylinder section is 
globally shrunk and a swelling appears in the region of 
the soliton. 
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